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Qh' Abstract 

"^ . The electronic structure, lattice dynamics, and electron-phonon coupling of the boron-doped 

diamond are investigated using the density functional supercell method. Our results indicate the 
1-^ boron-doped diamond is a phonon mediated superconductor, confirming previous theoretical con- 

Q ■ elusions deduced from the calculations employing the virtual crystal approximation. We show that 

o : 

the optical phonon modes involving B vibrations play an important role in the electron-phonon cou- 
pling. Different from previous theoretical results, our calculated electron-phonon coupling constant 



> 

■^ • is 0.39 and the estimated superconducting transition temperature Tc is 4.4 K for the boron doped 

V^ , diamond with 2.78% boron content using the Coulomb pseudopotential fi* = 0.10, in excellent 

o, 

TJ" ' agreement with the experimental result. 

o' 

C3 PACS numbers: 74.25.Kc, 74.25.Jb, 74.70.Ad, 74.62.Dh 
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Recently, Ekimov et al} reported the discovery of superconductivity in the boron-doped 
diamond synthesized at high pressure and high temperature. Their measurements showed 
that the boron-doped diamond with a hole carrier density of 5 x 10^^ cm~^ is a bulk, 
type-11 superconductor below the superconducting transition temperature Tc ~ 4 K. Using 
a simplified McMillan formula^ by setting the Coulomb pseudopotential /x* to zero, they 
estimated the electron-phonon coupling (EPC) constant A = 0.2, indicative of weak EPC 

Before that, superconductivity in some doped semiconductors has been discovered 
experimentally,—'^ and superconductivity in many-valley degenerate semiconductor has been 
also predicted theoretically.— As the first group- IV diamond type semiconductor with super- 
conducting properties, the boron-doped diamond superconductor has been studied by several 
groups. Baskaran et air- suggested that when the boron doping concentration increases to 
a critical value ric ~ 4 x 10^^/cm'^, an Anderson-Mott insulator to resonating valence bond 
(RVB)i superconductivity transition takes place. By first principles density functional per- 
turbation calculations!^ employing the virtual crystal approximation (VCA),^ Boeri et al^ 
substantiated that the recently discovered superconductivity below 4 K in 3% boron-doped 
diamond is caused by the coupling of a few holes at top of the cr-bonding valence band to the 
optical zone-center phonons, similar as in MgB2,'^^^^^*^^*^ albeit in 3 dimensions. Another 
first principles study employing the VGA and frozen phonon method also indicated that the 
EPC is the likely superconducting mechanism.— Regardless of their similar conclusion in 
these two studies, we note that there are some discrepancies in these two studies, especially 
in the EPC strength and T^: A = 0.27 and T^ = 0.2 K when /i* = 0.10 for 3% doping 
concentration in the first study, however, A = 0.55 and Tc = 9 K when /i* = 0.15 for 2.5% 
doping concentration in the second study. The difference in doping concentration makes the 
discrepancy even more prominent since larger doping concentration leads to larger A and Tc- 
Moreover, though qualitative agreement with experiment is found in these two studies, some 
obvious disagreements still exist, in particular, T^. in the first study is too small, and T^ in 
the second study is too large though one can obtain Tc = 4 K by setting /i* to a somewhat 
large value, i.e., 0.20. The disagreements between experiment and theory might indicate 
the inapplicability of the VCA to the newly found boron-doped diamond superconductor. 
In fact there are some cases where the VCA failed, e.g., the VCA was inappropriate for 
calculating the band gaps of GaPN and GaAsN,^^ and it described badly the structural 
and electronic properties of the quaternary alloy GaAlAsN.i^ On the other hand, though 



computational demanding, the supercell method is usually reliable.— li^ To bridge the gap 
between experiment and theory, and to verify the phonon mediated superconducting mech- 
anism in the boron-doped diamond, here we report a first principles supercell calculation 
on the EPC of the boron-doped diamond. Our results support the conventional phonon 
mediated superconducting mechanism in boron-doped diamond and the calculated Tc is in 
excellent agreement with the experimental value. 

Electronic structure calculations and geometrical optimizations are performed using 
density functional theory (DFT)^^*^^ in the local density approximation (LDA) .^°i^^ The 
electron-ion interaction is described by ultrasoft pseudopotentials,— which allow a low cutoff 
energy (25 Ry in this work) in the plane-wave expansion. The phonon and EPC calculations 
are carried out using density functional perturbation theory in the linear response.^ Within 
the phonon mediated theory of superconductivity, Tc can be estimated using McMillan's 
solution of the Ehashberg equation,— 

r-<exp/ i-m + ^) 1 (1, 

°"1.20 P\ A-f,-(l+0.62A)J' ^' 

where A is the EPC constant, uf^ is the logarithmically averaged characteristic phonon fre- 
quency, and /i* is the Coulomb pseudopotential which describes the effective electron-electron 
repulsion. The EPC constant A is calculated as an average over the A^ q— points mesh and 
over all the phonon modes, A = J^au '^<iu/N, where Aqy is the electron-phonon interaction 
for a phonon mode u with momentum q. The modes responsible for superconductivity can 
be identified from the Eliashberg function. 
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Then wj^ is calculated as. 



ujf^ = exp s -r dua'^F{u) \nuj/uj > . (3) 

We use the supercell technique to model the boron-doped diamond. In order to study the 
dependence of the EPC on the doping concentration, we choose two models: 2x2x2 and 
3x3x2 diamond supercells with a carbon atom substituted by a boron atom, named models 
I and II respectively. The total B content (Cb) for experimental samples is 2.8±0.5%, which 
is smaller than the total B content (6.25%) in model I and very close to the total B content 



(2.78%) in model II. The two models are shown in Fig. QJ The optimized lattice constant 
for diamond is 3.57 A, agreed well with the experimental lattice constant (3.566 A)i and 
previous theoretical results.— For the boron-doped diamond we check that our calculations 
reproduce the slight lattice expansion, less than 0.3%.^ Since the not very large boron doping 
has negligible effect on the lattice constant of diamond, we use our optimized lattice constant 
(3.57 A) in all subsequent calculations. 

The k-point integration for geometrical optimization, construction of the induced charge 
density, and calculation of the dynamical matrix is performed over a4x4x4 (2x2x3) 
Monkhorst-Pack grid^^ for model I (II), and a finer 8x8x8 (4x4x6) grid is used in 
the phonon linewidth calculations where the convergence in the k-point sampling is more 
difficult than that for the phonon calculations. The dynamical matrix and phonon linewidth 
are computed ona3x3x3 (2x2x3) q-point mesh for model I (II) , and a Fourier interpolation 
is used to obtain complete phonon dispersions. 

The Lowdin population analysis has been carried out to get the local density of states 
(LDOS). The total electronic density of states (TDOS) and LDOS for models I and II are 
plotted in Fig. El (The TDOS for the undoped diamond is also shown as a reference). The 
TDOS clearly indicates a degenerate or metallic behavior in both models, contrasting sharply 
to the semiconductor behavior of the undoped diamond. The TDOS and LDOS for model 
I are very similar with those for model II except that the width of the acceptor bands (the 
bands between the valence top and the fermi level) for model I is larger than that for model 
II. For both cases, the width of the acceptor bands is larger than that for a boron-doped 
diamond 64-atom supercell,— suggesting a dependence of the width of the acceptor bands 
on the boron doping concentration. When the doping concentration decreases, the width 
of the acceptor bands decreases. At very small doping concentration, the acceptor bands 
even no longer overlap with the valence band edge of the diamond resulting in a threefold 
degenerate acceptor state with a hole bind energy of Eb ~ 0.37 eV.-^ Since there are only s 
and p electrons in boron doped diamond, the electron correlation should not be very strong. 
The metallization in the boron-doped diamond is reasonablely caused by the increased boron 
content, doesn't necessarily resort to the Anderson impurity model , ^^i^^ which is suited for 
describing strong correlated systems. Fig. Efb) and (c) contain the LDOS plots of the B 
atom and the average LDOS of the C atoms. From Fig. |21we can clearly see that the LDOS 
around the Fermi level for B is larger than that the average LDOS for the C atoms. However, 



it doesn't mean the electronic states near the Fermi level are localized around B since the 
sum of the LDOS for all C atoms is larger than the LDOS for B. 

The calculated frequency of the highest optical phonon at T for the updoped diamond is 
1295 cm^^, agreed well with previous LDA calculations,— i^^ a little smaller than experimen- 
tal frequency, 1332 cm~^.^ The phonon band structures for the undoped and boron-doped 
2x2x2 diamond supercell are shown in Fig. Efa) and (b) respectively. Here, we plot 
the phonon band structure for the diamond supercell just for comparison with that for the 
boron-doped diamond. The phonon band structure for the diamond supercell is more com- 
plex than that for the unit cell since there are much more carbon atoms in the supercell. 
There are more branches for the boron-doped diamond since some phonon degenerates are 
removed due to the symmetry breaking, especially the highest optical phonon at F is a one- 
fold Ai mode, as shown in the inset of Fig. El^b). Beside that, a general effect to the phonon 
band structures resulting from the boron doping is the phonon softening. Especially the 
softening of the optical phonons is sizeable, e.g., the highest frequency of the zone center 
optical phonons decreases 99 cm~^ to 1196 cm~^ for model I. Comparing our results with 
those obtained from the VGA, we find that the softening of the highest optical phonon 
modes at F in our supercell calculations is much smaller (the decrease in theirs is 265 cm~^ 
for diamond with 5% boron conten1>i^). Zhang et al. observed that the zone-center optical 
phonon line at 1332 cm~^ downshifted for boron doped diamond film.^^ Part of the reason 
for the downshifting might be the softening of the optical phonon due to the EPC We can 
see that there is a slight upturn of the uppermost mode especially when moving from F to 
X for both updoped and doped diamond, similar with the previous results.—"^ Moreover, 
the upturn for the boron-doped diamond is more noticeable due to the larger EPC for the 
zone center optical phonons. 

The partial phonon density of states (DOS) for atom a is defined as: Paiyj) = 
Z]qSj=i |6a(q5 j)P'^(^ — "^(q^j))) where A^ is the total number of atoms, q is the phonon 
momentum, j labels the phonon branch, ea(q, j) is the phonon displacement vector for atom 
a, and u;(q, j) is the phonon frequency. The total phonon DOS, B partial phonon DOS and 
average partial phonon DOS of all C atoms for both models are shown in the upper panels 
of Fig. m The shapes of the phonon DOS for these two models are very similar, while the 
phonon softening effect in model I is stronger due to the larger boron concentration. The 
lower panels of Fig. EJare the Eliashberg function a^F{uj) for both models. The extremely 



weak signal in the low frequency part (lower than 80 meV) of a'^F{uj) indicates very weak 
electron acoustic phonon coupling. Detailed analysis shows that the B vibrations related 
EPC is the largest among all atoms in both models I and II (14.5% in model I and 9.3% in 
model II). Thus the B related phonon modes play an important role in the EPC due to the 
large B electronic LDOS near the Fermi level. As shown in the lower panels of Fig.^J a'^F{uj) 
has sizeable contributions from phonons with medium frequency in our results. However, 
Boeri et al. suggested that a^Fiuj) vanishes for phonon frequencies below that of the optical 
zone-center modes, then jumps to a maximum, and finally falls.— Their relative simple pic- 
ture for a^F{ijj) should stem from the fact that there are only three optical phonon modes 
in their virtual crystal calculations. Since we involve much more phonon modes to calculate 
a^F(co'), our results should be more reliable. 

The calculated TDOS at the Fermi level {N{Ep)), average frequencies uf^ , EPC con- 
stants A, and superconducting transition temperatures T^ for models I and II are listed in 
Table HI As expected, the N{Ef) for model I is larger than that for model II due to the 
larger boron concentration. The N{Ef) for the boron-doped diamond is consistent with 
the virtual crystal calculations: N{Ef) (in states/(spin ■ eV ■ diamond unit cell) is 0.062 
for 2.78% boron vs 0.060 for 2.5% boron,i^ 0.106 for 6.25% boron vs 0.08 for 5% boron.^ 
The smaller N{Ef) for the boron-doped diamond than MgB2 would result in weaker EPC 
in the boron-doped diamond. The large tof^ for the boron-dope diamond leads to high Tc, 
as demonstrated in McMillan's formula.— Thought the phonon softening increases with in- 
creased boron concentration, uf^ for model I is larger than that for model II, i.e., 1287 K 
vs 1218 K. The smaller tof^ for model II results from the substantial EPC in the medium 
frequency region. The EPC constant A for model I is 0.56, larger than that (A = 0.39) 
for model II. The typical value for /i* is in the range 0.10 to 0.15. The calculated Tc for 
model I is 23.6 (11.5) K for /i* =0.10 (0.15), and for model II it is 4.4 (0.9) K for /i* =0.10 
(0.15). The boron concentration for model II is very close to the experimental value and 
the calculated Tc agrees well with experimental Tc, about 4 K. We can see Tc increases with 
increased boron concentration due to the increased A. Detailed analysis reveals that the 
EPC is peaked at F, which is also found in previous study.— Lee et ali^ overestimated A 
and Tc since A should be averaged over the whole BZ, but their frozen phonon calculations 
only took into account the zone center phonons. 

Our treatment neglects the anharmonic corrections. For MgB2, Choi et alM showed that 



A is decreased from 1.0 to 0.78 after including partial anliarmonic corrections. However, 
more recently, Lazzeri et al^ suggested that the effects of anharmonicity on T^ and A in 
MgB2 are indeed negligible by explicitly taking into account the scattering between different 
phonon modes at different q points in the whole BZ. Boeri et al}^ considered the anharmonic 
corrections in the boron-doped diamond using frozen phonon calculations and found that 
the effect of anharmonicity is small, i.e., in general A decreases 0.03 after taking into account 
the anharmonicity. So our main results remain essentially unchanged even after including 
the anharmonic corrections since the anharmonic corrections for the boron-doped diamond 
are small, as discussed by Lee et al.;^ 

In summary, we have carried out a first principles study on the heavily boron-doped 
diamond. Optical phonons in diamond are softened after doping with boron. The boron 
related vibrational modes contribute an important part to the Eliashberg function a'^F{uj). 
Superconductivity in the boron-doped diamond is found to be mediated by the EPC. By 
using the supercell technique, we resolve the discrepancy between theoretical results based 
on the VGA and experimental data, and the calculated Tc is in excellent agreement with the 
experimental result. 
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TABLE I: Computed electronic DOS at the Fermi level {N{Ep)), average frequencies cof^^ , EPC 
constants A, and superconducting transition temperatures Tc for models I and II. The two values 
for Tc correspond to two different values of fi* (0.10 and 0.15). N{Ep) is in states/(spin • eV • 
diamond unit cell). Boron content Cb is also shown for each model. 



Cb JV(j^fX^(K) a T, (K) 

Model I 6.25% 0.106 1287 0.56 23.6,11.5 
Model II 2.78% 0.062 1218 0.39 4.4,0.9 




FIG. 1: (Color online) Structures of (a) model I and (b) model II. Refer to the text for the 
description for the models. 
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FIG. 2: (Color online) Electronic DOS for the undoped diamond and two boron-doped diamond 
models, (a) shows the TDOS of the undoped diamond, (b) and (c) are the TDOS and LDOS plots 
for models I and II. The LDOS for C shown here represents the average LDOS of all the C atoms. 
TDOS is in states/ (spin • eV • diamond unit cell). The unit of LDOS is states/ (eV • supercell). 
Energy is relative to the Fermi level Ep (for the undoped diamond, the valence top is taken as zero 
enegy point). 
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FIG. 3: (Color online) Phonon band structures for (a) 2x2x2 supercell diamond and (b) 2x2x2 
supercell diamond with a carbon atom substituted by a boron atom (model I). The insets show 
the high energy optical phonon branches near F both for the undoped and doped systems. 
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FIG. 4: (Color online) Total phonon DOS, boron partial phonon DOS and average partial phonon 
DOS for all carbon atoms, and Eliashberg function cP'Fiio) for (a) model I and (b) model II. The 
unit for the total phonon DOS is states/ (meV • diamond unit cell). Partial phonon DOS shown 
here is in arb. unit. 
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